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Abstract 

We show explicitly that the leading soft gluon pT distribution, predicted by 
Kovner, McLerran, and Weigert after solving classical Yang-Mills equations, 
can be understood in terms of conventional QCD perturbation theory. We also 
demonstrate that the key logarithm in their result represents the logarithm 
in DGLAP evolution equations. 

PACS numbers: 12.38.Bx, 12.38.Mh, 24.85. -hp 
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Typeset using REVT^ 



I. INTRODUCTION 



In ultra-relativistic heavy ion collisions, physical observables sensitive to a few GeV 
momentum scale, such as the mini-jet production, will be dominated by scattering of soft 
gluons from both heavy ion beams. Understanding the distribution of soft gluons formed 
in the initial stage of the collision is particularly interesting and important. In terms of 
conventional QCD perturbation theory, a calculable cross section in high energy hadronic 
collisions is factorized into a single collision between two partons multiplied by a probability 
to find these two partons of momentum fractions Xi and X2, respectively, from two incoming 
hadrons. The probability is then factorized into a product of two parton distributions 0(xi) 
and 0(x2), which are probabilities to find these two partons from the respective hadrons [|l|. 
Because of extremely large number of soft gluons in heavy ion beams, it is natural to go 
beyond the factorized single-scattering formalism to include any possible multiple scattering, 
and long range correlations between soft gluons from two incoming ions. 

Recently, McLerran and Venugopalan (MV) developed a new formalism for calculation 
of the soft gluon distribution for very large nuclei p|J^. In this approach, the valence quarks 
in the nuclei are treated as the classical source of the color charges. They argued that 
the valence quark recoil can be ignored in the limit when the gluons emitted are soft. In 
addition, because of the Lorentz contraction, the color charge of the valence quarks is treated 
approximately as an infinitely thin sheet of color charge along the light cone. With these 
assumptions, the gluon distribution function for very large nuclei may be obtained by solving 
the classical Yang-Mills Equation Using the classical glue field generated by a single 

nucleus obtained in the MV formalism as the basic input, Kovner, McLerran, and Weigert 
(KMW) computed the soft gluon production in a collision of two ultra-relativistic heavy 
nuclei by solving the classical Yang-Mills equations with the iteration to the second order 
0. The two nuclei are treated as the infinitely thin sheets of the classical color charges 
moving at the speed of light in the positive and the negative z directions, respectively. 
Following this approach, the distribution of soft gluons at the rapidity y and the transverse 
momentum px in nuclear collisions can be express as |^ 



where g is the strong coupling constant, A^c = 3 is the number of the color, and ^cutoff is 
a cutoff mass scale Note that there was a factor of vr misprint in Eq. (50) of Ref. [^, 
as pointed out in Ref. 0. In Eq. /i^ is the averaged color charge squared per unit area 
of the valence quark, and St is the transverse area of the nuclei. The /i^ and the St are 
related as 



where Ng is the number of valence quarks in the color charge source. The number distribution 
in Eq. can be also expressed in terms of the cross section 




(1) 
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In deriving Eq. the following relation was used 



da _ ^ dN 
dyd'^PT ^ dyd^pT 

The key result derived in Ref. Eq. (|l|) (or equivalently, Eq. (|D), is potentially very 
useful in estimating the production of mini-jet rates, and the formation of the possible 
quark-gluon plasma at RHIC . The purpose of this paper is to understand the respective 
role of perturbative and non-perturbative QCD in deriving the expression in Eq. (|I]) (or 
that in Eq. (^), and explore under what kind of approximation this result matches the 
conventional perturbative calculation. 

KMW's derivation for Eq. (0) is based on the following physical picture: in ultra- 
relativistic heavy ion collisions, gluons are produced by the fields of two strongly Lorentz 
contracted color charge sources, which are effectively equal to the valence quarks of two in- 
coming ions. In order to understand KMW's result in terms of the language of perturbative 
QCD, we consider a specific partonic subprocess: qq qqg, as sketched in Fig. |^. If we 
assume that the incoming quarks qq are the valence quarks in the initial color charge sources, 
the partonic subprocess in Fig. |1] mimics the physical picture adopted in KMW's derivation. 
However, as a Feynman diagram in QCD perturbation theory, the single diagram shown in 
Fig. |I| is not gauge invariant. As we demonstrate in Sec. under certain approximations, 
the contributions extracted from the diagram in Fig. || to the leading soft gluon production 
in Eq. is gauge invariant; and therefore, the physical picture proposed by KMW for soft 
gluon production is preserved. 

In Sec. m, within the framework of conventional perturbative QCD, we calculate the 
gluon production through the partonic subprocess qq —>■ qqg, as shown in Fig. 0, at the soft 
gluon limit. With our explicit calculation of this subprocess, we demonstrate that Eq. 
(or the cross section in Eq. (^) at A^g = 1 can be reproduced. Through our derivation, 
we show that the key logarithm in^p"^/ A^^^^jj-) in Eq. (|I]) (or in Eq. (§)) is basically the 
logarithm from the splitting of the incoming quark to the soft gluon in Fig. |I|. In terms of 
the conventional QCD factorization formalism such logarithm is normally factorized into 
the distributions of the coUinear gluons inside the incoming hadrons, and the logarithmic 
dependence of the distributions is a direct result of solving the DGLAP evolution equations 

i- 

Finally, in Sec. we discuss the relations between the MV formalism and that of 
the conventional QCD factorization. We explicitly demonstrate that KMW's result can 
be reduced to the factorized formula in the conventional perturbative QCD, if we replace 
the charge density for the classical color charge //^ (or equivalently Ng) by the valence 
quark distributions of the nuclei, and absorb the logarithm d-nip^ / Is?^^^^^ into one of the 
valence quark distributions. We point out that with the higher order of iteration, KMW's 
approach may include the multi-parton dynamics which is not apparent in the conventional 
perturbative calculation. 



II. FACTORIZATION AND GAUGE INVARIANCE 

As we discussed above, the partonic process qq —>■ qqg, as shown in Fig. |l|, mimics the 
physical picture adopted in KMW's derivation, if we assume that the incoming quarks qq 
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are the valence quarks in the initial color charge sources. However, in general, the Feynman 
diagram shown in Fig. |l| is not gauge invariant by itself. In this section, we discuss how to 
extract the gauge invariant leading contribution from the diagram in Fig. 1, and what kind 
of approximation we need to take in order to extract such leading contributions. 

For the production of gluons, we evaluate the invariant cross section, daqq^g/ dycfpr-i 
with y and px the repidity and the transverse momentum of the produced gluon. We label 
li and I2 as the momenta of the two incoming quarks, respectively, and we choose ki and /c2 
to be the momenta for the two gluons emitted from the initial quarks. We have = for 
the final-state gluon because of its on-shell condition. For the other two gluons, ki and ^2 
can not be on shell at the same time, because ki and /c2 come from different directions, and 

In general, the cross section can be written as 

X 

da = — IM]"^ dps , (5) 

where s = {li + and |Mp is matrix element square with the initial-spin averaged and 
the final-spin summed. In Eq. (H), dps is the phase space, and can be expressed as: 

dps = ^(2T)i((ii - - m') X ^(2,r)*(((, - k^f - m') 

where ^2 = p — ki, and m is the quark mass. For simplicity, we assume that both incoming 
quarks have the same mass. In high energy collisions, we set the mass of light quarks to be 
zero. Because of the gluon propagators, as shown in Fig. |l|, the matrix element square |Mp 
has the following pole structure: 

poles = , „ ^ . , „ ^ . , ^ ^ . , ^ ^ . . (7) 

kf + le kf — le K2 + le K2 — le 

When integrating over the phase space, we see that the leading contribution comes from 
the terms with k\ ^ Q 01 {p — kiY = k\ limit. As pointed out above, k\ and k^ can 
not be zero at the same time. Therefore, to calculate the leading contribution, we can first 
calculate the diagram in /cj ^ limit. The total leading contribution is just twice of it, 
because the diagram is symmetric for ki and ^2. 

When we take ki 0, the integration become divergent. Therefore, an introduction of 
a cutoff is necessary for obtaining a finite contribution from the diagram in Fig. |l|, and the 
corresponding contributions are sensitive to the cutoff. To derive the leading contribution 
at fc^ — > limit, we perform the collinear approximation ki ^ xli + 0(A;it), with kiT ~ 
^cutoff << Pt, where Acuta ff is a collinear cutoff scale 0. This approximation means that 
the leading contribution is from the phase space where almost all transverse momentum of 
the final-state gluon comes from the gluon of ^2, and ki is almost collinear to li. After 
such collinear approximation, the cross section in Eq. can be approximately written in 
a factorized form ilH : 
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where the overall factor of 2 is due to the fact that the leading contribution come from two 
regions corresponding to A;^ ^ and 0, respectively. In Eq. (||), Pi-^^ki{x, kiT < Pt) 

represents the probability of finding an almost collinear gluon with the momentum fraction 
X from an incoming quark of momentum /i, and 

Pz,^.,(x, k,T <Pt) = j ^,x5{x- |) |M,_,p (27r)5((/i - hf - m') . (9) 

The diagram for \Mg_^g\'^ can be represented by Fig. H{xli, h.p) in Eq. (|^) is effectively 
the hard scattering between the gluon of ki = xli and the incoming quark of l2- and given 
by 

H{xh,k,p) = H{xh,h,p) (27r)5((/2 + xh~pf - m^) , (10) 

where H{xli, h-iP) is given by the diagrams shown in Fig. ^ 

In addition to the diagram in Fig. |l|, in general, we also need to consider the radiation 
diagrams shown in Fig. ^. Similarly, the contribution of Fig. ^ and Fig. ^ can also be 
written in the same factorized form: 

with i = a,b. Here P/^^fc^(x, kiT < Pt) is defined by Eq. (||). Ha{xli, hjP) and Hb{xli, h^p) 
are the hard scattering parts from the diagrams in Fig. ^ and Fig. ^d, and they are repre- 
sented by Fig. ^ and Fig. |^, respectively. With the contribution from Fig. ^a and Fig. ^b, 
Eq. (ID changes to 

X [H{xli, l2,p) + Ha{xli, l2,p) + Hb{xli, hyP) + interference terms] 

+0(^^H^f^) , (12) 
Pt 

with the approximation ki = xli + 0(^1^). 

Feynman diagrams shown in Fig. ^ and Fig. |^ form a gauge invariant subset for calculat- 
ing the hard scattering parts, if(x/i, /2,p)'s in Eq. (|12|). Fig. |^ and Fig. Pd are effectively 
the s-channel and w-channel diagrams for the gq gq partonic process. Since we are only 
interested in the soft gluon limit, when |t| << s, the contribution from these two diagrams 
can be neglected, in comparison to the contribution from the diagram in Fig. 0. In addition, 
under the collinear expansion ki ^ xli, the gluon line which connects the partonic parts 
Pi^^ki and H{xli,l2,p) is effectively on the mass-shell, and therefore, the partonic parts, P 
and H in Eq. (^ are separately gauge invariant. 

Similar arguments can be held for the situation when fc^ ~ or k2T << Pt- For example, 
after the collinear expansion for k2, the contributions from diagrams shown in Fig. ^ and 
Fig. ^ can be neglected in the soft gluon limit. 
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Therefore, with the approximation of fc^ ~ (or ~ 0) and the soft gluon hmit, and 
a proper choice of the gauge, the dominant contribution for the partonic process qq 
qqg comes from the diagram shown in Fig. |1]. In the next section, we derive the leading 
contribution of the partonic process qq qqg with the above approximation. 



III. DERIVATIONS 

Following the discussion in last section, we now turn to explicit calculation of the leading 
contribution to the gluon production from the partonic process qq — >■ qqg, shown in Fig. |l|. 
As shown in Eq. (^), the leading contribution of this partonic process can be factorized into 
two parts: Pi^^^kiix, kiT < Pt) and H{xli, h^p)- Ph^ki{x, kir < Pr) represents the splitting 
of the quark to the soft gluon with momentum fraction x, and H{xli,l2,p) represents the 
scattering between the gluon of momentum xli and the other incoming quark of momentum 
I2. Pi-i^^ki{x, kiT < Pt) and H{xli, h-iP) are represented by the diagrams in Fig. ^and Fig. |^, 
respectively. In the following derivation, we choose the center of mass frame, with 

h = /i-, It) = (/+, 0, 0) and h = (0, /„, 0) . (13) 

The definitions of the plus and minus components of the four momentum p = {po,Pi,P2,P3) 
are: 

P0+P3 P0-P3 

We also introduce two useful vectors, n and n as: 

n = (0,l,Or), n = (l,0,OT). (15) 

As we discussed above, Pi-^^^kiix, kiT < Pt) and H{xli, h^p) are separately gauge invari- 
ant. To derive the complete leading contribution, we choose ■ A = gauge to calculate 
Pii_>fci(a;, kiT < Pt) and n ■ A = Q gauge for calculating the H{xli, h.p)- From Eq. (^ and 
the diagram shown in Fig. ^ we have in n ■ A = gauge, 

Pi^^k. (x, k,T < Pt) = C,^, g^l^^x6{x- ^){2n)6{{h - k,)' - m') 

x^Tr(7 ■ hn ■ ih - k,)^^)^^^^{-gn. (16) 

where Cq^g is the color factor. In Eq. (|^), the gluon polarization tensor is defined as 

P.,(fci) = -g^, + + (17) 

The four dimension integral d^ki = dki+dki^ ■ ndklj.. We can use the 5-function 6{x — y^) 
to fix ki+, and use 5((/i — ki)^ — w?) to fix ki_. We have 
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Substituting Eq. (|T8D into Eq. (0), and working out the trace, we obtain 



Pi^^kiix, kiT < Pt) = C, 



1 + (1 -X 



1 2 



X 



cutoff i-1 



Nl - 1 ( 1 + (1 - x) 



in 



X 



A2 



p±_\ 

cutoff 



(19) 



Choosing n ■ A = gauge, we derive the partonic scattering part of the H{xli,l2,p) 
defined in Eq. (10) from the diagram shown in Fig. |^, 

Hixh, h,p) =g'\ d^, Tr(7 ■ /27^'7 " (k - h)^') 



V2 rv2 

-2xli + p) + (-2p + x/i) V" +{p + xhfg"^ 



x[{p-2xhYg''^ +{xh+prg 



up 



-2p + xhYg 



u pal 



where d^iy is defined as 



dpu g^ii/ ~l~ ^fifT'u ~\~ ^fj.fT'u ■ 

In Eq. (pop, the gluon polarization tensors are given by 

PpAp, n) = -gp„ + — , 

p ■ n 

Ppp>[k2,n) = -gpp> + - 



k2 ■ n 

Using the relations k^ = {p — xliY and p^ = 2p^p_ — p"^ = 0, we have 

k'2 'p^p 

Substituting Eq. (p3D into Eq. (0), and after some algebra, we obtain 



I p± 




\xl+^ 


1 Pi-^ 



H{xhj2,p) = V (xs- 2xl+p_Y + i2xl+p_)i2p+p. 



(20) 



(21) 



(22) 



(23) 



(24) 



where s = (/i + I2Y = 2/+/_ is the total invariant mass squared of the two incoming quarks. 
Substituting Eq. (^^ into Eq. (|20|) , and after taking into account of the color factor 1/2, we 
obtain the hard scattering function H{xli,l2,p) as 



P± 

xlj 



Pt s — 2/+J9. 



-6{x 



Hixhj2,p) = (27r)V 

\ z / \ xi 

X (xs - 2x/+]9_)^ + (2x/+p_)(2p+p_) 
Combining Eq. (p^ , Eq. (|25|) and Eq. (H), we obtain 



2p+l- 



s — 2/+p_ 



(25) 
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(Pp (27r)4 \2N, 



X / ixSix - 1(1 + (1 - x]"- I ''^ 



2/+p_ J s{s — 2/+p_) 



X 



X 



We define the soft gluon hmit as 



and (27) 



At this soft gluon limit, we have s = 2/_|_/_ >> 2/+p_; and from the 5-function in Eq. (p6D, 
we have 

X f« < 1, and 1 + (1 - a;)2 ^ 2 . (28) 
Substituting Eq. ( pH]) into Eq. and taking the soft gluon limit, we obtain 



Defining y = ^in{^^), we can rewrite the cross section in terms of variable y 



This is our final result. Eq. ( pUD shows the same pj^ dependence as the result obtained by 
Kovner, McLerran, and Weigert 0. The difference between Eq. (^ is the factor of N^. 
Eq. (^) is obtained by calculating the leading contribution of the subprocess qq — > qqg, for 
which Nq effectively equals to one. If we consider the total number of the quarks in the 
charge sources of both sides, we need to multiply A^^^ to Eq. (|30|) , and our result reproduces 
the result obtained by KMW. 



IV. SUMMARY AND DISCUSSIONS 

In this section we discuss the similarities and differences between KMW's result, Eq. (|l]) 
or Eq. (^, which was obtained in McLerran- Venugopalan formalism, and our result, Eq. (pOD, 
which was obtained in the conventional perturbative QCD formalism at the leading loga- 
rithmic approximation. 

Our result can be reexpressed in terms of the usual factorized cross section in perturbative 
QCD. When we consider the collision between two nuclei, we can treat the two incoming 
quarks in Fig. |I] as coming from two nuclei respectively. In this picture, the number of the 
valence quark Nq is replaced by the quark distribution in the nuclei. The cross section in 
Eq. (ID (or equivalently Eq. (p9|)) is just the partonic cross section for the collision between 
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two quarks. In terms of the parton model, the cross section between the two nuclei A and 
B can be expressed in the following form: 



E 



dZi dZ2 fq/A{Zl) fq/B{Z2) E 



da, 



11^9 



d^p 



(31) 



Here zi and Z2 are the momentum fractions of the quarks, and fg/A^Zi) and fq/B{z2) are the 
quark distributions (or quark number densities) of the two nuclei. If we denote pa and as 
the momenta for the two nuclei respectively, then zi = h/pA and Z2 = h/vs- Substituting 
Eq. (D into Eq. (|l]), we have 



E 



d^p 



2(27r)3 2S 



f — — f — fq/A{zi)fg/B{z2)Ph^kAxi, kiT < Pt) H{xih, l2,p) 

fq/A{Zl)fq/BiZ2)Pl^^k2{X2, hr < Pt) H{li,X2l2,P) (32) 

Xo 



where the overall factor 2 in Eq. (|^) is now represented by the two terms, and S = {pa + 
PbY ~ 2pA ■ Pb- In Eq. (|32D, Xi = ki/k with i = 1,2, and h = zipA and I2 = Z2Pb- If we 
denote the momentum fraction of gluon ki with respect to pa as z'-^ = ki/pA, and ^2 with 
respect to ps as z'2 = ^2/^5, we can rewrite Eq. (|32D in terms of z'l and z'2: 



E- 



da 



AB^g 



d^p 



2(27r)3 25 

X 



dz[ dz'2 



Z^ Zo 



I J ^ fq/A{zi) Pi,^kAz'i/zi,kiT < Pt) 



X 



2(27r)3 25 



/ — fq/B{z2) P/2^fc2 (4/2:2, k2T < Pt) 
J Z2 

^) J~r^ / ^ /9M(2;i) Pi,^kAz'/zi, kiT < Pt) 



H{z[pa,z'^Pb,p) (33) 



X fq/Biz) H{z'pA,zpB,p) 



(34) 



In obtaining Eq. (p4D, we used the fact that the partonic scattering part H{ki,k2,p) in 
Eq. ( p3D is symmetric under the exchange of the ki and ^2 at the soft gluon limit. 

According to the QCD factorization theorem [Q], we see that the part inside the square 
brackets is actually the gluon distribution from nuclei A (or B) at the factorization scale 



= p^, with only the quark splitting function |T0 



— fq/A{zi)Ph^kAz'i/zi, kiT < Pt) 

Z\ 

+term from gluon splitting. 
Using Eq. (^5]), we can then reexpress Eq. (|33D as: 



(35) 
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,daAR^n 1 1 /■ dz' dz 

z 

+fq/A{z)fg/B{z',fJ,l = Pt)H{zpa,z'pb,p) , (36) 



which is the factorized formula for two-to-two subprocesses in the conventional perturbative 
QCD for the nucleus-nucleus collisions. In KMW formalism, only the valence quark color 
charge was used as the source of the classical charge of colors. As a result, the gluon splitting 
term in Eq. (|35| ) is neglected for the distribution fg/A- 

Our discussions above show that the soft gluon distribution in heavy ion collisions ob- 
tained in KMW's approach can be understood by calculating the partonic process qq qqg. 
To relate KMW's result to the factorized formula in the conventional perturbative QCD, we 
need to: (1) replace the charge density for the classical color charge /i^ (or equivalently Ng) 
by the valence quark distributions of the nuclei; (2) absorb the logarithm inlp'^/Al^^^j^j) 
into one of the valence quark distributions, which effectively becomes the gluon distribution 
of one of the initial nuclei. From Eq. (|33|) , we identify that the ln{p^ / K^^^^^j j) factor in 
Eq. (|I|) (or in Eq. (^) comes from the logarithm of the splitting of the incoming quark 
to the soft gluon in Fig. 0. In terms of the conventional QCD factorization theorem [|l|, 
such logarithm is normally factorized into the distributions of the collinear gluons inside 
the incoming hadrons, as demonstrated in Eq. (^), and the logarithmic dependence of the 
distributions is a direct result of solving the DGLAP evolution equations |Q. 

From the above comparison, we conclude that by solving the classical Yang-Mills Equa- 
tion to the second order in iteration, KMW's result reproduces the result of conventional 
perturbative QCD at the leading logarithmic approximation, with the convolution over the 
parton number densities inside the nuclei replaced by the effective numbers of the valence 
quarks. The logarithmic dependence shown in KMW's result basically describes the loga- 
rithmic DGLAP evolution of the quark distributions. However, in addition to the valence 
quarks, the glue field at small x can be produced by all flavor partons that have larger 
momenta. 

The McLerran-Venugopalan formalism was later further developed by Ayala, Jalilian- 
Marian, Kovner, McLerran, Leonidov, Venugopalan, and Weigert |[TT|Jl^ . The major im- 
provement to the McLerran-Venugopalan model is to include the harder gluons into the 
charge density /i^ and treat the charge source as an extended distribution which depends on 



the rapidity These improvements lead to the "renormalization" of the charge density. 



It was showed that the renormalization group equation for the charge density can be re- 



duced to the BFKL equations in some appropriate limits fl^ . It will be very interesting 
to see if KMW's approach, after including higher orders of iteration, can show the parton 
recombination ||T^ and other non-perturbative effects which are not apparent in the normal 
perturbative calculation. 
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FIG. 1. Square of the leading Feynman diagram to the process: qq qqg. 
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Fig. 2 

FIG. 2. Feynman diagram for the sphtting of g — > 5. 
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Fig. 3 

FIG. 3. Leading Feynman diagram contributing to the hard scattering part H{xli,l2,p)- 
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(a) (b) 
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(c) (d) 



Fig. 4 



FIG. 4. The rest Feynman diagrams to the process: qq qqg, in addition to the diagram 
I. 1. 




(a) (b) 
Fig. 5 

FIG. 5. Feynman diagrams contributing to Ha{xli,l2,p) (a), and Hi,(xli,l2,p) (b). 
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